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ABSTRACT 


This  technical  note  presents  an  improved  Kaiman  filter  excisor  for  dealing  with  liarrow- 
band  Gaussian  noise.  For  this  type  of  interference,  both  envelope  and  phas4?  vary  with 
time.  Since  the  Kalman  filter  approach  is  close  to  optimum  for  interferers  with  conslani 
or  very  slowly  varying  envelopes,  performance  for  the  narrowband  Gaussian  noise  case 
will  be  significantly  less  than  optimum.  Ways  and  means  of  improving  the  fK*rformance 
by  up  to  5  dB  are  presented.  These  include  the  use  of  an  alternate  filtering  approach  to 
estimate  the  envelope  to  counter  group  delay  effects,  and  phase-smoothing  to  suppress 
some  of  the  additional  noise  caused  by  temporary  loss  of  lock  by  the  Kalman  filter  during 
those  instants  when  the  envelope  changes  sign. 

RESUME 

Cette  note  techn  que  pr^nte  un  algorithme  base  sur  I’utilisation  d’un  filtre  de  Kalman 
pour  I’excision  de  bruits  gaussiens  a  bande  etroite,  lesquels  possedent  une  enveloppe  et 
t  une  phase  variant  rapidement  dans  le  temps.  Etant  donne  que  le  filtre  de  Kalman  cst  quaru 

optimal  pour  des  signaux  interferents  a  enveloppe  variant  tr^  lentement,  la  performance 
,  diminue  considerablement  lorsque  I’interference  est  un  bruit  gaussien  a  bande  etroite. 

Une  methode  pour  ameliorer  la  performance  de  5  dB  est  proposee.  Cette  methodc  inclut 
I’utilisation  d’un  autre  filtre  pour  estimer  I’enveloppe  du  signal  dans  le  but  de  compenser 
I’effet  du  delai  de  groupe,  ainsi  que  I’utilisation  d’un  filtre  pour  adoucir  les  brusques 
changements  de  phase  causes  par  les  pertes  de  synchronisation  du  filtre  de  Kalman, 
lesquels  se  produisent  lorsque  I’enveloppe  du  signal  change  de  signe. 
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EXECUTIVE  SUMMARY 


This  technical  note  presents  an  improved  Kalman  filtering  approach  that  is  used  for  fil¬ 
tering  narrowband  Gaussian  noise  interference  out  of  direct  sequence  spread  spectrum 
signals.  These  signals  are  used  extensively  in  military  communication  systems.  The 
technique  described  herein  applies  equally  to  both  Electronic  Supfmrt  Measures  (KS.V!) 
systems  and  direct  sequence  spread  spectrum  communication  sytems.  In  the  former  appli¬ 
cation,  the  ESM  system  may  be  attempting  to  intercept  the  spread  spectrum  signal,  but 
the  narrowband  interference  rnay  be  hampering  this  cffoit.  In  the  latter  appiicatiuii,  uie 
spread  spectrum  communication  system  may  require  additional  assistance  to  suppress  the 
interference.  Since  the  open  literature  has  been  devoted  to  this  latter  case,  the  material 
presented  here  focuses  on  this  application. 

One  of  the  attributes  of  direct  sequence  spread  spectrum  communication  systems 
is  their  ability  to  combat  interference  or  intentional  jamming  by  virtue  of  the  system's 
processing  gain  inherent  in  the  spreading  and  despreading  process.  The  interference  can 
be  attenuated  by  a  factor  up  to  this  processing  gain.  In  some  cases  the  gain  is  insufficient 
to  effectively  suppress  the  interferer,  leading  to  a  significant  degradation  in  system  per¬ 
formance  aa  manifested  by  a  sudden  increase  in  bit  error  rate.  If  the  ratio  of  interference 
bandwidth  to  spread  spectrum  bandwidth  is  small,  the  interference  cam  be  filtered  out  to 
enhaince  system  performance.  However,  this  is  at  the  expense  of  introducing  some  dis¬ 
tortion  onto  the  signad.  This  process  of  filtering  is  sometimes  referred  to  as  interference 
excision. 

The  Kalman  filtering  approach  is  based  on  the  digital  phase-locked  loop  Kalman 
filter  aund  is  close  to  optimum  for  demodulating  am  FM-type  of  interferer.  Because  the 
interference  is  assumed  to  be  much  stronger  tham  either  the  signal  or  noise,  the  Kalmam 
filter  locks  onto  the  interference  and  produces  an  estimate  of  the  phase  amd  envelope  of 
the  interference. 

In  the  case  of  narrowband  Gaussian  noise,  the  envelope  varies  significantly.  The 
time-vairying  nature  of  the  envelope  leads  first  to  an  unknown  group  delay  in  the  envelope 
estimate  if  an  IIR  filter  is  used  and,  second,  to  sudden  phase  chamges  when  the  envelope 
changes  sign;  for  this  latter  situation,  the  Kalmam  filter  loses  temporary  lock  in  its  phase- 
tramking  and  must  re-aicquire  the  interference.  Both  the  group  delay  and  phase- tracking 
problems  contribute  to  an  excess  amount  of  residuaJ  interference  at  the  output  of  the 
interference  canceller,  resulting  in  an  overadl  degradation  in  performance.  To  combat 
the  unknown  group  delay  caused  by  an  IIR  filter,  a  linear  phase  FIR  filter  with  known 
group  delay  is  used  in  the  envelope  estimator.  To  combat  the  excess  noise  caused  by 
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the  sudden  phase  changes,  additional  phase  smoothing  is  incorporated  in  tlie  Kalman 
filter.  Through  computer  simulation,  it  is  shown  that  both  enhancements  contribute  to 
an  overall  improvement  in  performance  of  4  to  5  dB  for  interference  bandwidlhs  ranging 
from  1%  to  5%  of  the  chip  rate. 
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1.0  INTRODUCTION 


Several  reports  have  been  written  by  the  author  on  the  subject  of  Interference  suppre.ss)<jn 
of  narrowband  interference  from  direct  sequence  spread  spectrum  communications  signals 
using  a  Kalman  filtering  technique  [1,  2,  3).  The  first  two  reports  dealt  with  inlerferers 
which  were  of  the  constant  envelope  type  for  which  the  Kalman  filtering  technique  provides 
close  to  optimum  performance.  The  last  report  evaluated  the  Kalman  filter  excisor  for  tin- 
case  of  narrowband  Gaussian  noise — a  more  difficult  type  of  interference  because  buih  ifse 
envelope  and  phase  vary  with  time.  This  technical  note  presents  an  improved  excisor  to 
deal  with  the  more  difficult  narrowband  Gaussian  noise  interference  problem  as  discussed 
in  [31- 

The  outline  of  the  technical  note  is  as  follows.  Section  2.0  presents  tlie  com¬ 
munications  model  used  in  the  simulations.  Section  3.0  describes  the  basic  interfered  e 
estimator  which  was  used  in  the  previous  reports  [1,  2,  3]  and  to  whicli  enhanceiiient.s  will 
be  made  to  improve  the  performance  for  narrowband  Gaussian  interference.  .Section  5.0 
presents  the  problem  to  be  solved.  These  include  two  things:  (a)  the  elfecls  of  group  delay 
caused  by  the  low  pass  filter  used  to  estimate  the  envelope;  and,  (b)  the  additional  phase 
noise  caused  by  a  loss  of  phase-lock  in  the  Kalman  filter  when  the  envelope  changes  sign. 
Enhancements  to  the  system  are  presented  and  evaluated.  Finally,  Sextion  5.0  couchnU's 
the  technical  note,  suggesting  areas  for  further  research. 

2.0  COMMUNICATIONS  MODEL 

The  basic  elements  of  the  BPSK  PN  spread  spectrum  receiving  system  are  shown  in 
Fig.  1.  The  received  waveform  r(t),  consisting  of  a  spread  spectrum  signal,  additive 
white  Gaussian  noise,  and  narrowband  interference  is  applied  to  a  bandp^iss  filter  with 
the  transfer  function  //(,p(/),  whose  output  is  defined  a.s 

u(0  =  5(0  +  n(0  +  i(f).  (1) 

The  bandpass  filter  for  the  application  considered  here,  is  assumed  to  be  a  filler 

matched  to  a  chip  and  centered  at  the  carrier  angular  frequency  u>o  of  the  spread  spectrum 
signal.  The  spread  spectrum  signal  is  defined  as 


s{t)  =  a{t)  cos{u)ot) 


(2) 
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Figure  1;  Spread  spectrum  communications  model 


where 

a{t)  =  '£D,b,{t-kn).  (:i) 

k 

In  Eq.  (3),  Dk  is  a  sequence  of  data  bits  of  amplitude  (±1)  and  duration  Tf,  seconds,  and 
bk(t  —  kTk)  is  the  PN  sequence  pattern  for  the  bit,  i.e., 

L 

^k{t)  =  ^Ckjq[t-  jTc)  (1) 

>=» 

with  L  being  the  number  of  pseudo  random  chips  per  bit,  or  the  processing  gain,  Ck^  is 
the  code  sequence  for  the  bit,  and  q{t)  represents  the  basic  chip  pulse  of  energy  E^. 

The  noise  n[t)  is  Gaussian  and  has  a  power  spectral  density 

Sn(/)  =  (5) 

where  NqI2  is  the  power  spectral  density  of  the  assumed  white  Gaussian  noise  from  the 
channel.  The  band  of  interference  is  defined  as 

i{t)  = /(t)cos(wof +  5(0)  (6) 

where  I{t)  is  the  interference  amplitude  and  5(0  is  the  phase  modulation.  It  has  been 
assumed  that  the  effect  of  the  bandpass  filter  Hbp{f)  is  negligible  on  the  interference  i(0. 

Referring  to  Fig.  1,  the  output  n(t)  of  the  bandpass  filter  Hbp{f)  is  bandpa.ss 
sampled  and  applied  to  a  limiter/bandpass  filter  and  interference  estimator. 

Consider  the  bandpass  sampler  first.  The  analog  signal  u{i)  from  Eq,  (1)  is 


2 


sampled  at  /,  =  2Rc  {rnf,  =  uJol'lir  +  Rcl2  for  some  integer  in)  wh<*re  Rc  is  the  chip  rate 
of  the  spread  spectrum  signal.  The  resultant  sampled  signal  is,  therefore, 

u„  =  3„  +  n„  +  i„,  (7,1 


where  s„  consists  of  the  sequence  {. . .  ,0,  (  — l)'*a„, 0, (  — 1  . .  .}  where  tiie  a,,  are 

of  energy  Ec  and  coded  according  to  CkjDk  for  the  chip  in  the  transmitted  bit. 
Hn  ‘  are  uncorrelated  Gaussian  noise  samples  of  variance  <7^^^  =  A'o/2),  and  in  is  liie 

sampled  version  of  Eq.  (6).  The  samples  u„  are  applied  to  the  interh’rerur*  eslimat(.»r  and 
interference  canceller. 

The  interference  estimator  produces  an  estimate  i„  of  the  interference  which  is 
removed  from  the  sampled  input  u„.  The  output  of  the  summer  is  decimated  and  sign 
reversed,  resulting  in  a  baseband  error  signal  e,.  This  error  signal  is  correlated  witli  ilie 
PN  sequence.  The  output  of  the  correlator  is  integrated  and  bit-detected. 

Consider  now  the  branch  containing  the  limiter.  Here,  is  appilied  to  a  llm- 
iter/bandp<uss  filter.  The  input  to  the  limiter  referenced  to  the  interference  is  reiiefined 
as 

“n  =  -/[fn  +  +  a',  H-  +  a'j  „]*  cospoii  -f  +  0„.nj  (8) 


where 


<Au,n  =  arctan 


^2.n/v^  + 


(2) 


+  ”i.n/V2  -f  a'j  „  j 

is  a  noise-like  phase  fluctuation  on  the  interferer’s  phase  and  is  due  to  the  noise  and 
spread  spectrum  signal.  The  terms  ^2.ni  ^i.n>  in-phase  and  quadrature 

components  of  the  noise  and  spread  spectrum  signal  with  respect  to  the  interference  phase 
On. 

The  output  of  the  limiter/bandpass  filter  is  [4] 


.  4/1'  ,  „  ,  , 

OPn  =  - cos  Ulon  -f  t/„  -4-  <Pu.n  , 

K 


(10) 


where  A'  is  the  limiter’s  output  level.  This  signal  is  redefined  as  (letting  .4'  =  v^^/4  for 
convenience) 

6p„  =  Spi  cos[won  +  +  4>u,n]-  (11) 

It  should  be  noted  that  for  large  interference-to-noise  ratios  in  which  the  interference  is 


'Coherent  bandpass  sampling  has  been  assumed,  so  that  the  in-phase  component  is 
(>^i,n/\/2)  co8(n7r/2)  and  the  quadrature  component  is  (nj.n/v^) 8in(n»^/2) 
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of  constant  envelope,  4>u,n  in  Eq,  (11)  is  approximately  Gaussian  [4].  The  siunphd  signal 
in  Eq.  (11)  is  what  is  processed  by  the  Kalman  filter,  which  estimates  the  j)hii.se  0,,  of  liie 
interference.  The  phase  estimate  is  denoted  as 

3.0  INTERFERENCE  ESTIMATOR 

The  interference  estimator  shown  in  Fig.  1  is  detailed  in  Fig.  2.  The  Kalnian  algoriliun 
[3]  produces  a  signal 

=  >/2cos(u;en  +  (12) 

which  can  be  used  as  a  basis  for  estimating  a  sampled  version  of  the  envelope  of  i(t)  li  e., 
/„  in  Eq.  (6)).  The  output  of  the  first  multiplier  is  a  baseband  term  and  is,  using  Fq.  (7j. 
excluding  the  y/2  factor, 

/'.,  =  InCOs{0n  -  0r,\n-l) 

+  n„  cos(won  +  )  +  a„  cos(t?„|„_i ) 

+  a„cos(2u;on  +  0„|„_i)  (13) 

Equation  (13)  consists  of  four  terms:  the  first  term  is  related  to  trie  desired  envclcjpe  of 


Figure  2;  Block  diagram  of  the  interference  estimator. 

the  interference;  the  second  is  approximately  Gaussian  baseband  noise  [5];  and  the  third 
and  fourth  terms  are  noise-like  terms  emanating  from  the  spread  spectrum  signal.  The 
fourth  term,  because  of  the  sampling  rate  conditions  discussed  in  [6],  is  essentially  filtered 
out  by  the  lew  pass  filter  of  bandwidth  Blpf  <  0.50  Hz  and,  thereforf',  will  be  ignored  in 
the  baseband  simulations  to  be  discussed  in  the  next  section.  The  term  /'/\/2  is  filtered, 
resulting  in  the  estimate  of  the  interference  envelope,  /„/\/2.  Combining  this  with  in 
Eq.  (12)  and  shown  in  Fig.  2,  yields  the  estimate  of  the  interference 

=  /„  cos(a;on  +  , ) ,  (11) 
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which  is  subtracted  from  u„  as  illustrated  itj  Fig.  1. 


4.0  NARROWBAND  GAUSSIAN  NOISE  INTERFERENCE 

The  interference  that  is  considered  is  narrowband  Gaussian  noise.  In  [1.  2],  the  envt-lojx* 
was  constant.  For  narrowband  Gaussian  noise,  the  envelope  varic*s  its  well  as  the  ph;o>e 
There  are  two  issues  to  be  concerned  with  here  that  are  related  to  the  time  varying 
envelope.  First,  the  low  pass  filter  used  in  t^stimatjng  the  envelope  /„  in  Fig.  2  will 
introduce  a  lag  in  /„  thereby  increasing  the  amount  of  residual  interference  at  the  output 
of  the  canceller.  Second,  a  time- varying  envelope  will  have  the  impact  of  introducing 
larger  phase  errors,  particularly  when  the  envelope  gets  very  small,  since  the  Kaiinan 
filter  will  not  be  able  to  track  the  fast  phase  transitions  during  these  instances.  These 
effects  will  now  be  dealt  with  at  length. 

4.1  FILTER  LAG  EFFECTS 

First  consider  the  effect  of  the  lag  introduced  by  the  Butterworih  filter  of  Intndwidtli 
Blpf  in  estimating  In  in  Fig.  2.  The  effect  can  be  obstrrved  if  one  graph-s  the  degree 
of  interference  suppression  at  the  output  of  the  excisor  as  a  function  of  d,  tlu*  frequency 
deviation  constant  of  the  Kalman  filter,  for  a  particular  interference  bandwidth  JI,  for 
several  values  of  Dipp.  Interference  suppression  is  defined  in  (1]  as 

5  =  I01og(^/F.)  (15) 

where  Ain  =  *n  ~  in  is  the  residual  interference,  and  P.  is  the  power  in  the  interference. 
As  an  example,  Fig.  3  illustrates  the  degree  of  interference  suppression  as  a  function  of 
d,  when  the  interference  bandwidth  is  B,  —  0.02  Hz.  The  interference  was  generated  by 
applying  white  Gaussian  noise  to  a  fourth  order  lowpass  Butterworih  filter  resulting  in 
the  baseband  interference  i„  =  /„cos(0„).  Observe  that  as  BiPf  increases,  the  degree 
of  suppression  improves.  This  phenomenon  is  due  to  the  reduction  in  the  lag  as  Bipp 
increases.  It  is  more  clearly  illustrated  in  Figs.  4  and  6,  in  which  sections  of  the  envelope 
In  and  its  estimate  /„  have  teen  plotted  for  Blpf  =  0.20  and  0.40  Hz  for  the  case  when 
d  —  0.30  rad. /sec/volt  which  is  approximately  the  minimum  in  Fig.  3.  A  comparison  of 
Fig.  6  with  Fig.  4  reveals  a  reduction  in  the  delay  but,  as  to  be  expected,  an  increase 
in  envelope  noise.  The  effect  of  the  delay  is  further  illustrated  in  F'igs.  5  and  7  for  the 

A 

overall  interference  i„  and  i„  for  the  two  cases.  In  Fig.  5,  one  can  also  see  the  effect  of 
the  phase-noise,  since  the  envelope  noise  for  this  case  wiis  small  as  indicated  in  Fig.  4. 
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From  these  results,  and  the  results  of  (l,  2)  concerning  constant  envelope  interfer¬ 
ence,  there  are  conflicting  requirements.  The  results  from  the  constant  envelope  analysis 
suggested  the  use  of  as  small  a  bandwidth  for  the  lowpass  filter  as  possible,  whereas  here, 
as  large  a  value  as  possible  is  required  to  reduce  the  lag,  but  at  the  expense  of  introducing 
more  envelope  noise  and  signal  distortion.  This  suggests  that  one  could  compensate  for 
the  delay  or  use  a  linear  phase  FIR  filter  with  some  predefined  delay.  Either  of  these 
approaches  will,  of  course,  translate  into  the  requirement  for  a  similar  delay  in  the  PN 
correlator.  The  remainder  of  this  technical  note  will  consider  an  FIR  low  pass  filter  with 
pre-defined  delay. 

To  accomodate  this  change,  the  interference  esimator  of  Fig.  2  will  have  to  be 
modified  somewhat  to  account  for  this  type  of  filter.  The  modified  version  of  Fig.  2  is 
illustrated  in  Fig.  8.  The  FIR  low  pass  filter  has  a  delay  D,  which  requires  that  j/n  in 
Fig.  2  and  Eq.  12  must  be  delayed  by  Z)  as  well  to  yield 

yn-D  =  \/2co3[a.\)(n  -  D)  +  K-d]  (16) 

before  combining  it  with  the  envelope  estimate,  /„_p,  to  produce  the  interference  estimate 
in-D-  This  requires  too  that  the  input  signal  u„  =  Sn  -f  n„  +  in  be  delayed  by  D  samples 
before  removal  of  the  interference  estimate.  The  output  of  the  canceller  is 

Cn-D  =  «n-D  “  ^n-D 

—  -Sn-D  +  +  in-£>  -  in-D-  (17) 

This  error  signal  is  then  applied  to  the  PN  correlator  whose  PN  sequence  is  delayed  by  D 
as  well. 

A  comparison  of  the  suppression  performance  of  the  systems  in  Figs.  2  and  8  is 
now  presented. 

The  bandwidth  of  the  FIR  and  Butterworth  filters  for  this  case  was  Blpf  =  0.20 
Hz.  The  FIR  filter  had  7  taps  with  rectangular  weighting,  the  interference  bandwidth  was 
Bi  =  0.02  Hz  and  Ei/No  =  12  dB.  The  results  are  illustrated  in  Fig.  9  and  indicate  that 
the  FIR  filter  improves  the  amount  of  interference  suppression  by  2.7  dB.  Furthermore,  a 
comparison  of  Fig.  3  with  Fig.  9  shows  that  the  FIR  filter  of  bandwidth  0.20  Hz  provides 
an  additional  1  dB  of  interference  suppression  over  the  Butterworth  filter  of  bandwidth 
0.40  Hz.  This  gives  one  some  idea  of  the  degradation  in  performance  produced  by  the 
additive  envelope  noise  using  the  interference  estimator  in  Fig.  2  for  the  larger  bandwidth 
situation. 

Examples  of  the  envelope  and  interference  estimates  are  also  illustrated  in  Figs.  10 
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and  11.  Notice  in  Fig.  11  between  200  and  300  seconds  the  additional  noise  due  to  phase 
noise.  Finally,  Fig.  12  shows  the  suppression  level  for  several  interference  bandwidlhs, 
with  Blpf  =  0.20  Hz  and  EbfNo  =  12  dB.  It  is  important  to  note  here  that  as  the  band¬ 
width  of  the  interference  increases,  the  vadue  of  d  at  which  the  minimum  occurs  increases, 
emphasizing  the  requirement  for  a  larger  Kalmaui  filter  bandwidth,  which  is  controlled  by 
d  (everything  else  being  constant),  to  cope  with  the  faster  rate  of  change  of  the  phase. 

4.2  PHASE-NOISE  EFFECTS 

The  second  issue  to  consider  is  the  tendency  of  the  Kalman  filter  to  temporarily  lose 
track  when  the  phase  varies  rapidly  during  those  instances  when  the  envelope  becomes 
small.  It  is  conceivable,  therefore,  that  with  some  aulditional  phase-smoothing  apart  from 
that  provided  by  the  Kalman  filter,  the  phase  estimates  can  improve,  thereby  reducing 
envelope  and  phitse-noise  caused  by  inaccurate  phase  estimates. 

To  accomodate  this  idea,  an  additional  modification  must  be  made  to  the  inter¬ 
ference  estimator  in  Fig.  8;  this  change  is  illustrated  in  Fig.  13.  In  this  figure,  the  phase 
estimate  from  the  Kalman  filter  in  Fig.  4  of  [3]  is  applied  to  a  simple  .V-point 

smoother  implemented  as  an  FIR  comb  filter  that  introduces  the  same  delay  D  as  the 
FIR  filter  used  in  filtering  the  envelope.  The  result  is  the  smoothed  phase  which  i.s 
applied  to  a  phase  modulator,  yielding  the  signal 

yh~D  =  \/2cos(u;o(«  -  D)  +  K-d\-  ( 

This  signal  is  combined  with  the  envelope  estimate  /„_£>  to  produce  the  interference 
estimate 

in-D  =  L-d  cos[u;o(n  -  D)  +  (19) 

that  is  then  subtracted  from  the  delayed  input  u„_£>  to  yield  the  error  signal  Cn-D- 

Figure  14  shows  a  section  of  the  phase  Juid  its  Kalman  filtered  estimate 
for  the  case  when  the  interference  bandwidth  is  Bi  =  0.02  Hz  and  d  =  0.30  rad. /sec. /volt 
using  the  system  in  Fig.  8.  This  should  be  compared  with  Fig.  15,  which  is  a  7-point 
smoothed  version  of  <9n|n-i  based  on  the  system  in  Fig.  13;  the  effect  of  smoothing  the 
phase  estimate  is  quite  apparent.  Figure  16  shows  a  magnified  version  of  a  segment  of 
interference  i„  along  with  the  estimates  obtained  with  and  without  phase  smoothing.  Here 
too,  the  improvement  is  quite  noticeable. 

The  level  of  interference  suppression  as  a  function  of  d  is  illustrated  in  Fig.  17  for 
several  interference  bandwidths.  As  can  be  seen,  the  degree  of  suppression  has  improved 
by  3  to  4  dB  in  comparison  with  the  results  in  Fig.  12  in  which  no  phase-smoothing 
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was  used.  In  addition,  there  is  a  broader  range  of  values  for  d  over  which  the  maximum 
suppression  level  can  be  maintained. 

The  bit  error  rate  performance  of  the  two  interference  estimators  in  Figs.  8  and 
13  is  presented  next  bls  determined  from  computer  simulation.  The  test  conditions  were 
as  follows.  A  total  of  50,000  bits  were  transmitted  with  a  processing  gain  of  Z.  =  20, 
implying  a  total  of  1,000,000  chips.  The  interference  was  narrowband  Gaussian  noise 
ranging  in  bandwidth  from  5,  =  0.01  Hz  to  0.05  Hz.  The  approximated  optimum  values 
for  d  obtained  from  Figs.  12  and  17  were  used  in  the  Kalman  filter.  The  term  q/  was  set 
at  0.0003125  Hz  and  the  spectral  density  of  the  observation  noise  No(,,/2  was  obtained 
from  Eq.  (23)  of  [1],  with  replaced  by  /^,  the  mean-squared  value  of  the  envelope.  The 
bandwidth,  Blpf,  selected  for  the  linear  phase  FIR  filter  was  0.20  Hz  and  introduced  a 
delay  of  Z)  =  3  samples.  The  phase  estimate  ^n|n-i  in  Fig.  13  was  applied  to  a  7-point 
rectangular  smoother,  which  also  introduced  a  delay  of  D  =  3  samples. 


Figure  3:  Level  of  interference  suppression  as  a  function  of  d  for  several  lowpass  filter 
bandwidths.  The  interference  is  narrowband  Gaussian  noise  of  bandwidth  0.02  Hz. 
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The  results  of  the  two  interference  estimators  are  plotted  in  Figs.  18  and  19.  The 
bit  error  rate  performance  of  the  interference  estimator  with  phase-smoothing  (Fig.  19) 
indicates  a  significant  improvement  in  performance  over  its  non-smoolhed  counterpart 
shown  in  Fig.  18.  In  fact,  the  improvement  is  4  to  5  dB  in  each  of  the  cases  when  bit 
error  rates  at  Ei,INo  =  12  dB  are  compared.  These  results  also  indicate  the  role  of  the 
phase-noise  in  so  far  as  its  effect  on  the  bit  error  rate  for  narrowband  Gaussian  noise  is 
concerned.  This  was  also  exhibited  in  the  amount  of  interference  suppression  between 
the  two  approaches  (Figs.  12  and  17).  There  too  the  improvement  in  suppression  was 
approximately  4  dB. 


Figure  4:  An  example  of  a  portion  of  the  envelope  /„  and  its  estimate  /„  when  the  filter 
is  of  the  Butterworth  type  of  bandwidth  B^pf  =  0.20  Hz  with  d  =  0.30  rad. /sec. /volt 
amd  interference  bandwidth  Bi  =  0.02  Hz. 
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Figure  5;  An  example  of  a  portion  of  the  interference  in  and  its  estimate  when  the 
filter  is  of  the  Butterworth  type  of  bandwidth  Bipf  =  0.20  Hz,  d  =  0.30  rad. /sec. /volt 
and  interference  bandwidth  J5,  =  0.02  Hz. 
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Figure  6:  An  example  of  a  portion  of  the  envelope  /„  and  its  estimate  /„  when  the  filter 
is  of  the  Butterworth  type  of  bandwidth  Bipf  —  0.40  Hz  with  d  =  0.30  rad. /sec. /volt 
and  interference  bandwidth  Bi  =  0.02  Hz. 
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Figure  7:  An  example  of  a  portion  of  the  interference  i„  and  its  estimate  when  the 
filter  is  of  the  Butterworth  type  of  bandwidth  B^pf  =  0.40  Hz,  d  =  0.30  rad. /sec. /volt 
and  interference  bandwidth  Bi  =  0.02  Hz. 
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d  (rad./sec./volt) 


Figure  9:  Comparison  of  the  interference  suppression  levels  S  produced  by  a  7-tap  FIR 
filter  and  4**  order  Butterworth  filter  when  Bi  —  0.02  Hz,  Bipp  =  0.20  Hz,  and  Eb/No  = 


Time  (sec.) 


Figure  10:  An  example  of  a  portion  of  the  envelope  and  its  estimate  /„_o  when 
the  filter  is  of  the  linear  phcise  FIR  type  of  bandwidth  Bi,pf  =  0.20  Hz  with  d  =  0.30 
rad./sec./volt  and  interference  bandwidth  Bi  =  0.02  Hz. 
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Figure  11:  An  example  of  a  portion  of  the  interference  and  its  estimate  in-o  when 
the  filter  is  of  the  FIR  type  of  bandwidth  Blpf  =  0.20  Hz,  d  =  0.30  rad.  /sec. /volt  and 
interference  bandwidth  Bi  =  0.02  Hz. 
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Figure  12:  Interference  suppression  level  S  as  a.  function  of  d  for  several  interfere 
bandwidths,  Bi,  ranging  from  0.01  Hz  to  0.04  Hz,  with  Bipp  =  0.20  Hz. 
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Figure  15:  An  example  of  the  phase  0„  and  a  7*point  smoothed  version  of  its 

Kalman  filter  estimate  0n|n-i,  with  interference  bandwidth  Bi  —  0.02  Hz  and  d  =  0.30 
rad. /sec. /volt. 
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Figure  16:  An  example  of  a  segment  of  interference  with  and  without  phcise  smoothing 
interference  bandwidth  =  0.02  Hz  and  d  =  0.30  rad. /sec. /volt. 
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Ef,lNo  (dB) 


Figure  19:  Bit  error  rate  for  the  case  of  narrowband  Gaussian  noise  with  bandwidths  Bi 
ranging  from  0.01  Hz  to  0.05  Hz,  using  the  optimum  values  of  d  obtained  from  Fig.  17, 
and  an  FIR  filter  with  7-point  averaging  on  the  phase  estimate 


5.0  CONCLUSIONS 


This  technical  note  has  presented  an  improved  Kalman  filtering  excisor  to  deal  willi  the 
more  general  type  of  interference — narrowband  Gaussian  noise,  which  would  be  represen¬ 
tative  of  the  case  when,  say,  several  independent,  narrowband  interferers  are  lumped  close 
together  in  some  region  of  the  frequency  band  also  occupied  by  the  direct  sequence  spread 
spectrum  signal.  Under  these  conditions,  both  the  envelope  and  phase  of  the  composite 
interference  will  vary,  resulting  in  poorer  performance  of  the  basic  Kalman  filter  excisor 
described  in  earlier  reports. 

This  technical  note  presented  the  essence  of  the  problem  in  so  far  as  the  Kalman 
filter  is  concerned  when  narrowband  Gaussian  interference  is  present.  Two  enhancements 
were  suggested  and  evaluated.  It  was  shown  that  up  to  5  dB  of  improvement  could  be 
attained  with  the  revised  excisor. 

It  should  be  noted  that  only  one  value  of  low  pass  filter  bandwidth  D^pf  (i.e., 
0.20  Hz)  was  selected  for  these  comparisons,  and  only  one  type  of  phase-smoother.  Other 
combinations  may  produce  superior  performance;  however,  further  work  would  have  to 
be  conducted  to  assess  if  indeed  this  would  be  possible.  Finally,  given  the  significance 
of  the  phase-noise  at  the  output  of  the  Kalman  filter  for  narrowband  Gaussian  noise 
interference,  other  architectures  perhaps  can  be  developed  using  the  phase-smoothing 
concept  to  further  improve  performance. 
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